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ABSTRACT 


An  approximate  method  for  the  determination  of  tunnel  wall  interference  on 
two-dimensional  non-lifting  supercritical  airfoils  is  presented.  The  method 
assumes  classical  linearized  subsonic  interference  theory  appl i d to  an  air- 
foil whose  thickness  is  increased  so  as  to  match  the  far  field  transonic  vel- 
ocity potential.  The  interference  potential  at  the  wing  trailing  edge  is 
taken  as  the  parameter  controlling  the  shock  location  on  the  airfoil.  Sample 
computations  of  circular  arc  airfoils  in  tunnels  with  solid,  porous  and 
slotted  walls  are  shown  to  adequately  predict  the  shock  displacement  produced 
by  the  wall  constraints.  Also  presented  is  a numerical  analysis,  based  on  a 
shock  fitting  technique,  for  the  determination  of  exact  shape  and  location  of 
imbedded  shocks.  With  the  analysis  and  accompanying  numerical  code  some 
results  arc  provided  to  show  the  sensitivity  of  back  pressure  perturbations  on 
shock  location. 
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SECTION  I 
INTRODUCTION 


Since  the  early  I970's,  with  the  development  of  the  numerical  computation  tech- 
niques initiated  by  Murman  and  Coles\  it  has  been  possible  to  study  numerically 
the  transonic  flow  field  around  more  and  more  complex  configurations  in  free- 
air  as  well  as  within  the  confinements  of  wind  tunnels.  At  the  same  time,  a 
method  for  assessing  the  presence  and  level  of  wall  interference  in  transonic 

tests,  and  for  eliminating  this  interference  by  an  adjustment  of  wall  conditions, 

2 3 k 

was  suggested  by  Fern  and  Baronti  and  by  Sears  . More  recently,  Kemp  has 
introduced  the  concept  of  "correctable"  interference,  for  the  direct  determin- 
ation of  transonic  tunnel  corrections.  Thus,  by  a combination  of  analyses  and 
testing  techniques,  it  would  appear  that  significant  advances  can  be  expected 
in  developing  means  for  recognizing,  measuring  and  possibly  eliminating  tran- 
sonic tunnel  interference. 

Although  some  progress  in  wall  correction  techniques  has  been  reported,  notice- 
ably in  the  experiments  performed  at  Calspan^,  the  procedure  based  on  "adaptive" 
or  "self-correcting"  walls  appears  rather  complex  and  difficult  to  implement. 

On  the  other  hand,  systematic  numerical  computations  of  flow  fields  around  aero- 
dynamic configurations  \n  free-air  and  in  tunnels  can  indicate  the  presence 
of  interference.  The  computations  are,  however,  lengthy,  and  do  not  directly 
yield  the  level  of  v/all  interference  and  tunnel  corrections  as  in  classical 
subsonic  interference  theory.  The  direct  determination  of  wall  irterference 
and  tunnel  corrections  for  transonic  conditions  (subcritical  as  well  as  super- 
critical) v/ould  be,  on  the  other  hand,  very  useful  as  recognized  and  pursued 

678 

analytically  in  the  work  performed  by  Lo  and  Kraft  ’ ’ at  AEDC,  and  in  the 
work  performed  by  Kemp  , which  combines  numerical  computations  with  wind  tunnel 
pressure  measurements  at  the  wall  and  at  the  model  to  infer  tunnel  corrections. 

An  approach  which  contains  some  of  the  concepts  developed  by  Kraft  and  Lo,  and 
which  intends  to  provide  a direct  determination  of  tunnel  interference  by  a 
combination  of  approximate  analytical  methods  and  tunnel  measurements,  has  been 
the  subject  of  an  investigation  at  GASL  and  is  reported  here.  The  method  is 
based  on  the  observation  that  for  certain  wing  profiles  such  as  circular  arc 
airfoils  at  zero  angle  of  attack,  the  wall  interference  is  mostly  manifested  by 
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a shift  of  the  shock  position,  with  minor  alteration  of  the  wing  pressure  distri- 
bution upstream  and  downstream  of  the  shock,  and  by  a corresponding  variation, 
with  respect  to  free-air  conditions,  of  the  potential  at  the  wing  trailing  edge. 

This  variation  can  be  calculated  with  reasonable  accuracy  by  a standard  sub- 
sonic interference  method,  based  on  the  Prandt 1 -G1 aucrt  scaling  rule  and  applied 
to  a thicker  wing  profile  which  satisfies  the  far  field  transonic  velocity  poten- 
tial. The  method  is  applied  to  circular  arc  profiles  at  zero  angle  of  attack 
in  tunnels  v/i th  solid,  porous  and  slotted  walls  and  is  shown  to  provide,  in  a 
first  approximation,  a good  indication  of  wall  effect  on  shock  position. 

In  a parallel  effort,  also  described  in  this  report,  an  attempt  has  been  made 
to  assess  the  effects  of  back  pressure  on  pressure  distribution  and  shock  loca- 
tion on  an  airfoil.  The  effects  of  back  pressure  and/or  axial  pressure  grad- 
ients on  model  characteristics,  in  particular  on  shock  location,  have  been  rec- 
ognized repeatedly,  for  instance  in  connection  with  the  boattail  drag  of  axi- 
symmetric  models.  Similarly,  back  pressure  effects  were  held  responsbile  for 
the  displacement  of  shock  location  on  circular  arc  airfoils  from  expected  free 
flight  conditions,  in  a series  of  experiments  conducted  in  a joint  effort 
betv/ecn  GASL  and  WPAFB,'^’^'.  The  determination  of  the  effects  of  small  back 
pressure  perturbations  on  shock  location  docs,  however,  require  the  availabil- 
ity of  an  accurate  procedure  for  locating  the  shock.  In  this  respect,  an 
analysis  (and  a special  computer  code)  which  employs  a shock  fitting  technique 
and  which  satisfies  exactly  the  transonic  Ranki ne-Hugon iot  shock  Jump  lelation, 
has  been  developed  to  determine  the  exact  shape  and  location  of  an  embedded 

I 

shock.  j 

i 

The  analysis  is  an  extension  of  concepts  previously  developed  by  Murman  and  j 

by  Hafez  and  Cheng'^.  More  refinement  is  however  obtained  with  the  present  i 

procedure.  V/i th  the  shock  point  operator  of  Murman  the  shock  is  smeared  over 

several  grid  points;  with  the  shock  fitting  technique  of  Hafez  and  Cheng  the 

shock  jump  relations  are  applied  between  two  contiguous  points,  and  the  exact 

shock  location  is  not  provided.  The  present  code,  however,  determines  the 

exact  position  of  the  shock  betv/een  neighboring  grid  points  and,  thus  allows 

for  a more  precise  evaluation  of  the  effects  of  small  external  perturbation  on 

shock  location,  even  in  the  case  of  relatively  course  grid  networks.  A discus- 
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Sion  of  the  analysis  and  of  the  numerical  algorithm,  and  sample  computations 
showing  the  effect  of  back  pressure  on  shock  location,  are  presented  in  a 
following  section. 
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SECTION  II 


APPROXIMATE  ANALYSIS  OF  WALL  INTERFERENCE 


‘ EFFECTS  ON  SHOCK  LOCATION 

I 


I 

I Basic  Cons i derat  ions 

The  proposed  procedure  for  inferring  the  effects  of  wall  interference  on 
shock  location  is  based  on  a number  of  sequential  steps  each  of  which  may  con- 
tain certain  physical  assumptions  and  mathematical  approximations.  Each  step 
is  now  discussed  in  a sequential  manner. 

1.  A basic  assumption  underlying  the  procedure  is  that,  in  certain 
transonic  tests  of  non-lifting,  two-dimensional  airfoils,  wall  interference  will 
mainly  change  the  position  of  the  shock  and  only  sligiitly  the  distribution 

of  pressure  upstream  and  downstream  of  the  shock.  Such  a situation  is  typified 
by  the  result  of  Figure  (1),  taken  from  the  experiments  of  Reference  (Hi),  which 
show  the  pressure  distribution  over  an  airfoil  in  a tunnel  with  slotted  walls. 
The  data  along  curve  (a)  are  considered  interference-free  (small  model  with 
respect  to  tunnel  dimensions),  whereas  the  data  of  curve  (b)  are  affected  by 

; wall  interference  with  a consequent  upstream  shift  of  the  shock. 

2.  The  second  st(?p  of  the  procedure  is  to  postulate  that  the  shock 
position  on  the  airfoil  is  known  and  that  the  pressure  distribution,  before  and 
aft  the  shock,  can  be  cither  computed  by  an  approximate  analytical  method 

(the  method  of  Spreiter  and  Alksne'^  for  both  the  accelerating  and  decelerating 
portions  of  the  flov/  has  been  used  in  the  computations  performed  here  for 
circular  arc  airfoils),  or  that  this  pressure  distribution  is  the  one  directly 
measured  in  an  experiment.  (The  fact  that  experimental  results  may  contain 
wall  interference  does  not  detract  from  the  overall  logic  of  the  procedure.) 

Then,  by  integration  of  the  pressure  over  the  airfoil  for  several  shock 
locations,  a unique  functional  relationship  betv/een  the  shock  location,  x^ , 
and  the  difference  of  the  perturbation  velocity  potential,  between 

leading  and  trailing  edge  can  be  established  for  any  given  Mach  number  and 

j airfoil  geometry.  This  relationship  is  held  to  be  true  for  both  unconfined 

1: 

i and  confined  flows. 


I ^Further  qualifications  to  this  statement  will  be  provided  below. 
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figure  1.  PRESSURE  DISTRIBUTION  OVER  A 6%  CIRCULAR  ARC  AIRFOIL  IN  A 
SLOTTED  TUNNEL,  EXPERIMENTAL  DATA  FROM  REFERENCE  (]k) . 
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3.  The  next  step  is  to  calculate  the  transonic  far  field  potential  from 
the  known  pressure  distribution  over  the  airfoil.  As  has  been  established  in 
Reference  (1),  the  far  field  potential  of  a non-lifting  airfoil  is 


•{-(x.y) 


1 

271(5 


2 

X +3  y 


(1) 


where  the  doublet  strength 


D = 


CO 


F(fJdC  + ^ 


u^d^dC 


(2) 


I 

i 


1 


consists  of  a terni  proportional  to  the  airfoil  volume  and  of  a double  integral 
containing  the  nonlinear  contribution  of  the  flow  field  kinetic  energy. 


The  potential  given  by  Equation  (1)  has  the  standard  form  corresponding  to  a 
doublet  with  Prand 1 1 -Gi aue r t scaling.  Thus,  once  the  value  of  the  double  into- 

JL 

gral  is  determined,  a corresponding  body,  having  a "similar"  shape  F (fj  , can 
be  defined  to  yield  the  same  doublet  strength 


D = 2 


\ 


F*(C)dC 


(3) 


The  new  body,  which  evidently  has  a greater  volume  than  the  original  one,  will 
be  subsequently  used  to  determine  wall  Interference.  Before  doing  so,  a 
representation  of  the  velocity  in  the  flow  field  around  the  airfoil  must  be 
provided  in  order  to  compute  the  double  integral. 

k.  Since  the  early,  approximate,  analytical  methods  based  on  the  Green 
theorem,  such  as  the  method  of  developed  by  Spreiter  and  Alksne'^,  rational 
representations  of  the  mid  field  velocity,  u,  have  been  suggested.  Among 
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these  is  the  approximation  used  by  Kraft 


7 . 

1 n 


h i s 


integral  method, i.e., 


u(x,y)  = u(x,0)  e , 


(M 


relating  the  mid  field  velocity  to  the  velocity  on  the  airfoil,  with  the 
parameter  r(x)  satisfying  the  i rrotat ional i ty  condition  on  the  airfoil  and 
given  by 


r(x)  = - 

Ou/8y)^^Q 


(5) 


The  approximat ion  implies  that  u is  equal  to  7ero  in  the  regions  upstream  of 
the  leading  edge  and  downstream  of  the  trailing  edge.  Also,  it  cannot  be 
correct  for  negative  r,  since  it  v^ould  imply  an  expotential  increase  of  the 
perturbation  velocity  from  the  airfoil  outwards;  however,  by  taking  in 
Equation  (^)  the  absolute  value  of  r,  for  regions  where  r<0,  the  proper  decay 
can  be  recovered.  These  inaccuracies  are  recognized  as  conceptual  limitations 
ofthe  approximation.  The  errors  that  application  of  Equation  [^)  may  produce 
in  the  present  context  is  however  thought  to  be  acceptable  and  such  as  to  not 
invalidate  the  solution  procedure. 

5.  With  the  new  airfoil  of  shape  F (^)  , satisfying  Equation  (3),  sub- 
sonic linear  interference,  with  Prandtl-Glauert  scaling,  is  now  invoked. 

The  nev/  (thicker)  airfoil  is  assumed  to  be  placed  in  a wind  tunnel  with  either 
solid  or  ventilated  walls  (porous  and  slotted  of  infinite  length),  and  the 
perturbation  produced  by  the  wall  constraints  is  calculated. 

(Details  of  the  calculation  for  non-lifting  airfoils  are  given  in  the  Appendix.) 
This  value  of  is  then  added  linear!  ly  to  the  vabe  Aifi  determineed  in  Step  1, 
and  from  the  relationship  relating  to  A({),  the  new  position  of  the  shock 

is  determined. 

The  procedure  can  be  utilized  in  a direct  manner  to  determine  the  expected 
shift  in  shock  position  when  placing  an  airfoil  in  a wind  tunnel.  More  rewarding. 
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however,  is  the  application  of  the  procedure  to  wind  tunnel  airfoil  data. 
Inoeel,  by  the  experimentally  determined  pressure  distribution  on  the  airfoil 
it  will  be  possible,  by  following  the  steps  indicated  above,  to  determine  the 
value  of  Ai?'  induced  by  the  tunnel  wall  constraints  and  then  the  amount  of 
shock  displacement  caused  by  the  wall  interference. 


Application  to  a Circular  Arc  Airfoil 

The  pressure  distribution  over  a circular  arc  airfoil  under  supercritical  flow 
conditions  is  obtained  by  applying  the  local  linearization  method  of  Spreiter 
and  Alksne'^.  For  an  assumed  shock  position,  the  pressure  distribution  in  front 
and  behind  the  shock  is  obtained  by  applying  the  accelerating  and  decelerating 
solution,  respectively.  Accelerating  and  decelerating  solutions  are  matched 
at  the  shock  by  assuming  that  the  shock  is  normal  and  by  satisfying  the  normal 
shock  jump  relations.  Shov-/n  in  Figure  (2)  is  a comparison  between  a numerical 
solution  and  the  Sprei ter-Al ksne  solution  (with  the  shock  position  taken  to 
be  the  one  obtained  with  the  numerical  computation).  The  Spre i te r-Al ksne  solu- 
tion is  then  used  to  determine, in  function  of  free-stream  Mach  number,  the 
difference  in  potential,  A<}i,  between  leading  and  trailing  edge  in  terms  of 

shock  locations  x . Typical  results,  for  representative  free-stream  Mach  numbers, 

s • 

arc  plotted  in  Figure  (3). 


The  next  step  of  the  procedure  is  to  compute  the  far  field  potential.  Equations 
(1  and  2).  '‘or  a circular  arc  airfoil  of  thickness  x,  the  mid  field  velocity. 
Equations  (4  and  5),  can  be  v/ritten  as 


u = u 


where  u = As  shown  in  Figure  (4, a),  this  approximation  reproduces  very 

accurately  the  numerical  results  over  the  portion  of  the  wing  v;ith  negative  c^. 
In  the  regions  near  the  airfoil  nose  and  behind  the  shock,  whe'e  c^  is  positive, 
the  approximation  is  still  acceptable,  as  shown  in  Figure  (4,  ' , when  the  abso- 
lute value  of  Uq  is  taken.  The  strength  of  the  doublet  D is  then  provided  by 


^Numerical  solutions  have  been  obtained  with  the  GASL  transonic  program,  which 
includes  the  shock  fitting  technique  discussed  in  Section  III. 
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Spreiter  and  Alksne 


FIGURE  2.  PRESSURE  DISTRIBUTION  OVER  A CIRCULAR  ARC  AIRFOIL 
COMPARISON  BETWEEN  A NUMERICAL  SOLUTION  AND  THE 
SOLUTION  OF  SPREITER  AND  ALKSNE,  REFERENCE  (15). 
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the  relationship 


D = 2 


2t{1/A-x^)  dx  + (y+1 ) ! | e dydx 


2T  Y+1 
3 8t 


The  evaluation  of  the  integral  is  performed  numerically  and  the  value  of  D 
determined.  The  strength  of  tlie  doublet  is  nov-;  taken  to  correspond  to  that  of 
a new  circular  arc  airfoil  of  thickness  t for  which 

o-f' 


Hence,  the  new  airfoil  thickness  is  provided  by  the  expression 

. 1/2 


2 8t 


3 r 

2 8.  J 


The  airfoil  of  thickness  x satisfies  the  far  field  potential.  It  is,  however, 
of  interest  to  compare  with  the  numerical  results  the  normal  velocity  pertur- 
bation induced  at  closer  distances  from  the  airfoil  (in  correspondence  to  possible 

A 

locations  of  tunnel  walls)  by  an  airfoil  of  thickness  t . Such  a comparison 

is  shov/n  in  Figure  (5)  for  a circular  arc  airfoil  which  has  an  original  thickness 

of  6%  and  which,  with  the  inclusion  of  the  nonlinear  effects  embodied  by  the 

A 

double  integral,  has  a thickness  t =0.13-  The  analytical  results  with 

J. 

T =0.13  are,  necessarily,  antisymmetric  with  respect  to  the  airfoil  center  . 

(x=0  position)  in  contrast  to  the  actual  physical  situation  reflected  by  the 
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numerical  results.  The  level  of  the  perturbat ion  is,  hov/ever,  in  reasonable 
agreement  with  the  numerical  results  and  much  closer  to  these  than  the  analytical 
result  based  on  the  value  of  t=0.06. 

It  remains  now  to  determine  the  interference  potential  produced  by  the  ivall 
constraints.  Several  computations  have  been  performed  parametrically  in 
terms  of  free-stream  Mach  number  and  tunnel  height  for  tunnels  v/ith  solid, 
porous  and  slotted  walls.  For  a porous  or  solid  wall,  the  interference  poten- 
tial is  provided  by  Equation  (A, 8)  of  the  Appendix;  for  a slotted  wall  by 
Equation  (A, 13).  For  circular  arc  airfoils  the  integrals  of  Equations  (A, 9), 

(A, 10)  and  (A,1A)  can  be  calculated  in  closed  form.  The  integrals  (in  the 
variable  p)  in  (A, 8)  and  (A, 13)  have  been,  on  the  other  hand,  computed  numer- 
ically. Once  the  difference  of  the  interference  potential  betv/een  airfoil 
trailing  edge  and  leading  edge  is  computed,  i.e.  ' [_£  differ- 

ence is  added  to  the  value  of  pertaining  to  the  original  shock  location  x^ , 
and  the  nev/  shock  location  is  then  determined  from  Figure  (3)- 

Shown  in  Figure  (6)  is  the  effect  of  a solid  wall  on  shock  location  for  a 6^ 
thick  circular  arc  airfoil  at  M =0.91  in  a tunnel  of  height  twice  the  airfoil 
chord,  i.e.,  h/c=2.  (The  free  air  results,  taken  as  a reference,  are  those 
given  by  Kraft^.)  Also  shown  in  the  figure,  for  comparison,  are  the  results 
obtained  by  Kraft^  by  the  integral  equation  method  based  on  the  solution  of 
the  transonic  small  disturbance  equation.  Reasonable  agreement  between  the 
present  results  and  those  of  Kraft  can  be  observed.  Shown  in  Figure  (7)  are 
typical  results  pertaining  to  the  same  airfoil  in  a tunnel  with  porous  walls 
for  the  two  values  of  T=l/R=l<  and  10.  Good  agreement  with  the  results  of 
Kraft'  can  again  be  observed.  Finally,  the  results  given  in  Figure  (8)  refer 
to  the  same  airfoil  in  a slotted  tunnel  with  h/c=l  .5  and  2.00  and  with  a slot 
parameter  K=0.5.  As  indicated  in  Reference  (I4)i  a value  of  K=0.5  corresponds 
to  the  slotted  tunnel  of  the  Flight  Dynamics  Laboratory,  Wr ight-Patterson  air 
Force  Base,  the  facility  from  which  the  experimental  data  presented  previously 
in  Figure  (1)  were  taken.  Thus,  the  results  of  Figure  (8),  when  referred  to 
the  experimental  data  of  Figure  (1),  seem  to  provide  again  a reasonable  pre- 
diction of  shock  displacement  (upstream  displacement  in  this  particular  case) 
under  the  effect  of  wall  interference. 
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The  validity  and  extension  of  the  method  to  other  configurations  depends  evident- 
ly on  the  validity  and  generality  of  the  assumptions  that  have  been  introduced. 
One  assumption  is  that  wall  interference  does  not  nuch  change  the  pressure 
di  s t r i but  i 0[i  upstream  of  the  shock.  The  results  of  Figure  (9),  taken  from 
Reference  (7),  which  shov/  typically  the  effects  of  a solid  wall  on  the  pressure 
distribution  over  a NACA  12  airfoil,  seem  to  indicate  that  tlie  assumption  has 
not  a general  validity  and  that  the  wall  interference  on  the  pressure  distri- 
bution is  configuration  dependent.  However,  even  under  tlie  conditions  of  the 
figure,  it  seems  that  the  major  effect  of  wall  interference  is  the  shift  of 
the  shock  and,  thus,  that  the  procedure  can  be  perhaps  invoked  for  an  assess- 
ment of  first  order  corrections  for  general  non-lifting  airfoils. 


/ 
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FIGURE  9.  EFFECT  OF  A CLOSED  WIND  TUNNEL  ON  PRESSURE 

DISTRIBUTION  OVER  A NACA  0012  AIRFOIL,  REFERENCE  (7) 
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SECTION  111 

SHOCK  FITTING  TECHNIQUE  AND  DETERMINATION  OF  BACK 
PRESSURE  EFFECTS  ON  SHOCK  LOCATION 


Background 

During  the  past  several  years,  difference  formulas  have  been  successfully  em- 
ployed to  numerically  approximate  the  differential  equations  governing  a tran- 
sonic flow  field.  Successive  line  overrelaxation  (SLOR)  is  a reliable  and 
efficient  procedure  to  solve  the  resulting  set  of  simultaneous,  algebraic 
equations,  and  computer  codes  using  such  relaxation  schemes  have  been  devel- 
oped to  treat  the  transonic  small  disturbance  equation.  If  the  flow  is  super- 
critical these  codes  make  no  distinction  for  grid  points  in  the  vicinity  of  the 
shock  but  rather  allow  the  shock  wave  to  naturally  develop  in  the  course  of  the 
iteration  procedure  as  a discontinuity  smeared  over  a few  mesh  points.  Results 
with  embedded  shocks  do,  however,  indicate  a shock  pressure  rise  on  the  surface 
consistently  less  than  the  theoretical  value. 

1 2 

Murman  attempted  to  resolve  the  discrepancy  by  the  introduction  of  a shock 
point  operator.  At  a supersonic-subsonic  transition,  the  x difference  operator 
at  the  first  subsonic  point  dov-vns t rcani  of  the  shock  is  replaced  by  the  sum  of 
the  x-di f fercnccs  for  Ihe  elliptic  and  hyperbolic  operators.  Murman  then 
demonstrated  that  the  correct  normal  shock  jump  condition  occurs  over  three 
mesh  intervals.  For  the  case  of  an  infinite  oblique  shock,  the  shock  point 
operator  can  reproduce  the  correct  shock  jump  condition  at  upstream  and  down- 
stream flow  positions  which  are  several  grid  points  removed  from  the  shock 
position  (under  the  stipulation  that  the  oblique  shock  passes  through  two  grid 
points,  i.e.,  (dx/dy^=&x/NAy , where  N is  an  integer,  or  Ax, Ay  arbitrary  in  the 
limiting  case  of  an  infinitesimally  refined  grid  network). 

The  shortcomings  of  Murman's  shock  point  operator  are  discussed  by  Hafez  and 
Cheng'^,  viz.  (1)  the  requirement  that  grid  arrangement  be  such  that 
Ax/Ay=N(dx/dy)^  is  not  readily  obtainable  for  an  arbitrary  grid,  (2)  the  dis- 
continuity in  4i  across  the  shock  is  computed  in  a manner  that  will,  in  most 
cases,  underestimate  ^ on  the  upstream  side  of  the  shock  and  (3)  the  correct 
supersonic-subsonic  transition  cannot  be  accomplished  in  a single  jump  but 
requires  a few  grid  points. 
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Hafez  and  Cftcng'^  suggested  an  improved  shock  fitting  technique,  which  satisfies 
the  Rank i ne-[|ugon iot  condition  and  which  defines  more  exactly  the  shock  position 
even  with  coarser  grids.  The  shock  fitting  procedure  modifies  the  difference 
equations  mainly  in  the  treatment  of  grid  points  near  the  shock.  The  original 
difference  equation  is  replaced  by  the  difference  form  of  the  jump  relations 
at  the  point  on  each  horizontal  line  nearest  to  the  shock  on  the  downstream 
side,  and  backward  extrapolation  is  used  to  evaluate  the  potential  at  the  near- 
est grid  point  along  each  vertical  on  the  upstream  side.  For  a plane-normal 
shock,  the  above  shock  fitting  procedure  collapses  to  Murman's  shock  operator 
formulation.  In  general,  however,  the  two  procedures  do  not  agree  because  of 
the  difference  in  which  the  jump  in  transverse  velocity  i{>  is  computed. 


The  shock  fitting  procedure  of  Hafez  and  Cheng  does,  however,  apply  to  the  shock 
jump  conditions  between  neighboring  grid  points  rather  than  at  the  actual  shock 
position.  It  is  the  present  intent  to  further  refine  the  shock  fitting  tech- 
nique, removing  this  discrepancy  by  calculating  the  shock  position  within  the 
iteration  procedure  as  a discrete  discontinuity  across  which  the  transonic 
approxi ma t i on  of  the  Rankine-Hugoniot  conditions  apply. 

• Formulation 

The  transonic  small  disturbance  approximations  of  the  Rankine-Hugoniot  condi- 
tions are 

<4' 

■ (I'M  ) + ('i'  +'}'  ) = ■ - y (6, a) 

X 


> 

X 


dy  c 


(6,b) 


where  < > signifies  the  jump,  and  the  subscripts  1 and  2 represent  the  upstream 
and  downstream  side  of  the  shock,  respectively.  The  first  equation  is  the  tran- 
sonic shock  polar  and  the  second  equation  expresses  the  continuity  of  tangential 
velocity  across  the  shock.  The  latter  condition  can  be  immediately  satisfied 
by  requiring  the  potential  to  be  continuous  across  the  shock.  The  line  relax- 
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ation  procedure  to  v/h i eh  the  shock  fitting  method  will  be  applied  is  that 

11,12 


of  Murman 
that  of  Hafez  and  Cheng 


The  derivation  of  the  shock  fitting  procedure  closely  resembles 

13 


For  a supersonic-subsonic  shock  transition,  the  modifications  in  the  differ- 
ence equations  are  found  mainly  in  the  treatment  of  the  grid  points  on  each 
horizontal  line  nearest  to  the  shock  on  the  downstream  side,  and  of  the  near- 
est grid  point  along  each  vertical  across  the  shock.  These  points,  labelled 
P and  Q,  respectively,  are  shown  in  Figure  (10).  Three  types  of  local  shock 
inclination  indicated  in  Figure  (10)  nujst  be  distinguished  in  treating  point  P. 
By  examining  the  elliptic/hyperbolic  nature  of  the  grid  points  about  P the 
type  of  local  shock  inclination  is  determined.  The  shock  inclination  having 
been  classified  as  type  (a),  (b)  or  (c) , its  position,  6 priori  known  to  lie 
within  the  interval  P P,  is  calculated  by  the  following  procedure.  The 
interval  P P is  subdivided  into  tl-1  equally  spaced  intervals  as  shown  typically 
in  (b)  of  Figure  (10)  where  P , n=l , 2,  3 • • • , N is  a general  point 

I n 

vyithin  P P.  Assuming  the  point  P to  correspond  to  the  shock  position,  values 

n 

of  the  potential  and  of  the  x-velocity  component  at  the  upstream  side  of  the 
shock,  (<(>|)p  ^'^iV  ’ ' y ’ obtained  by  forward  extrapola- 

tion from  upstream  gri;i  points.  Continuity  of  the  potential  across  the  shock 
yields  the  potential  downstream  side  of  the  shock,  whereas 

Equation  (6, a)  yields  the  ^x-vclocity  component,  (o2)p  on  the  downstream 


side  of  the  shock 
side  of  the  shock 


The  first  ordei  Taylor  expansion  about  the  downstream 


♦p  = -('(Xp  ,yp)  + (<>^)p  (Xp-Xp  ) 

n n n 


(7) 


supplies  a potential  at  P,  denoted  by  The  error  function 

E(xp  ) = Jp  ■ 4’p^^  where  4’p^'^  is  the  potential  at  P at  the  previous 
n 

I 

iteration,  is  formed  and  the  position,  Xp  =x^ , in  the  interval  P P,  where 
the  error  function  is  zero,  is  identi f ied^as  the  shock  location. 


To  compute  the  RHS  of  Equation  (6, a),  values  obtained  at  the  previous  iteration 
step  are  used,  with  the  if  's  and  the  's  evaluated  by  central  differences  at 
A and  B,  and  D and  C,  respectively,  in  ^ the  manner  shown  in  (a),  (b)  and  (c) 
of  Figure  (10) . 
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The  difference  scheme  at  P is  written  so  as  not  to  cross  the  shock  wlicn 
approximating  derivatives  by  finite  differences.  The  original  difference 
equation  for  the  upstream  point  Q must  also  be  redefined  in  cases  (a)  and  (b) , 
because  one  of  the  five  points  in  the  backv/ard  hyperbolic  operator  lies  on 
the  other  side  of  the  shock.  The  needed  data  arc  supplied  by  forward 
(dov/nstream)  extrapolation  to  P from  three  upstream  points. 

For  a supercritical  flow  the  elliptic  and  hyperbolic  regions,  identified  by 

|_  “ “>  xj>0elliptic, 

are  separated  by  a sonic  boundary  and  by  a shock.  For  flows  accelerating 
through  a sonic  boundary  the  difference  form  of  the  equations  of  the  original 
relaxation  scheme  is  retained.  For  flows  decelerating  through  a shock,  however, 
the  difference  equations  of  the  original  scheme  are  replaced  by  the  difference 
formulation  described  above.  Following  Reference  (13).  to  distinguish  the 
shock  from  the  sonic  line,  ip^  is  required  to  decrease  by  a prescribed 
fraction  of  the  grid  size  in  a supersonic-subsonic  transition. 


Numerical  Results 

The  present  shock  fitting  procedure  has  been  investigated  intensively  and 
numerical  results ^have  been  compared  with  those  obtained  by  applying  the  method 
of  Hafez  and  Cheng. (A  numerical  code  based  on  the  analysis  of  Hafez  and  Cheng 
has  also  been  constructed  to  allow  direct  comparison  of  the  two  methods.)  A 
typical  comparison  is  provided  by  Figure  (11)  which  shows  the  pressure 
distribution  over  a 6Z  thick  circular  arc  airfoil  at  = 0.87*  The  grid 
network  consists  of  Al  x 20  points.  The  tv/o  sets  of  results  are  those  obtained 
after  1200  relaxation  steps.  Shown  in  Figure  (12)  is  the  algorithm  providing 
the  shock  location  within  two  neighboring  points.  Plotted  in  the  figure  is 
the  error  function  E(C^)for  two  points  above  a 6%  thick  circular  arc  airfoil 


The  potentials  at  P , P and  S are  used  to  form  the  difference  approximations 
of  th«'  x-de  r 1 va  t i ves , whereas  central  differences  at  C (for  shock  inclination 
types  (a)  and  (b)  of  '"igjre  (10))  or  at  P (for  shock  inclination  type  (c)  of 
Figure  (10)  ) are  used  to  approximate  the  y-derivati ves. 
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FIGURE  11.  COMPARISON  BETVEEN  PRESENT  SHOCK  FITTING  SOLUTION 
WITH  SOLUTION  BY  HAFEZ  AND  CHENG,  REFERENCE  (13) 
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FIGURE  12.  DETERMINATION  OE  SHOCK  POSITION  BETWEEN 
ADJACENT  GRID  POINTS 
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at  M = 0.94.  The  shock  location  is  identified  by  the  values  of  C for  which 

s 

E{C^)  = 0.  (?  is  the  stretched  x-coordinate  in  the  computational  space.) 

Systematic  calculations  have  then  been  performed  with  the  program  to  ascertain 
the  effects  of  back  pressure  on  airfoil  characteristics  and,  in  particular, 
on  shock  location.  The  representative  calculations  reported  here  consider 
a 6%  thick  circular  arc  airfoil  at  M = 0.9'  and  assume  several  constant  values 
(different  for  the  free  flight  values)  of  the  pressure  along  the  normal 
direction  at  3 station  approximately  one  chord  downstream  of  the  airfoil. 
Unperturbed  potential,  iji  = 0,  has  been  assumed  to  prevail  at  upstream  infinity 
and  at  an  infinite  distance  above  the  airfoil.  Clearly  indicated  in  Figure 
(13)  ai'c  a forward  shift  of  the  shock  for  values  of  back  pressure  larger  than 
the  free-air  values  and  a shift  toward  the  wing  trailing  edge  for  values  of 
C smaller  than  free-air  values.  The  effect  of  back  pressure  is  also  felt 
in  the  region  above  the  airfoil  as  shown  in  Figures  (lA)  and  (15),  which 
present  pressure  distributions  and  flow  deflection  angles  along  a line  at  a 
distance  of  y/c  = I.58  above  the  airfoil.  However,  the  pressure  distri- 
butions on  the  airfoil  upstream  of  the  shock  is  only  slightly  perturbed  with 
respect  to  the  free-air  value. 

• 

As  indicated  by  the  result  of  Figure  (l3)i  the  shift  in  shock  location  with 
a perturbation  of  back  pressure,  provides  the  major  contribution  to  the 
variation  of  the  potential  at  the  wing  trailing  edge.  The  availability  of 
approximate  analyses,  such  as  the  one  described  in  the  previous  section,  to 
account  for  back  pressure  effects  (or  similar  perturbations  from  free-air 
values  dov/nstream  of  a model)  by  a corresponding  variation  of  trailing  edge 
potential  would  then,  again,  be  very  useful. 


FIGURE  15.  FLOW  DEFLECTIONS  ALONG  LINE  y/c  - 1.58  FROM  AIRFOIL  FOR  SEVERAL  VALUES  OF  BACK  PRESSURE 
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In  transonic  testing  it  would  be  useful  to  possess  the  capability  of  correcting 
the  tunnel  data  in  the  san'e  fashion  followed  in  subsonic  testing.  Attempts  are 
now  being  made  to  develop  capabilities  that  permit  one  to  quantify  transonic 
wall  interference;  the  concept  of  correctable  interference,  which  combines 
tunnel  measurements  with  numerical  analyses,  is  an  attempt  in  this  direction. 

The  method  suggested  here  is  an  even  simpler  procedure  for  an  approximate  deter- 
mination of  tunnel  interference.  Based  on  the  observation  that  a n 'or  mani- 
festation of  wall  interference  is  a shift  of  shock  position,  a methodology  is 
presented  for  the  determination  of  the  shift  of  the  shock  due  to  wall  constraints. 

The  feasibility  and  relative  accuracy  of  the  method  is  proved  herein  only  for 
simple  configurations  such  as  a circular  arc  airfoil  at  zero  angle  of  attack. 

An  extension  of  the  concept  and  of  the  method  to  other  non-lifting  configura- 
tions such  as  airfoils  v/i th  blunted  leading  edges,  and,  possibly,  to  lifting 
airfoils,  could  constitute  a worthwhile  effort.  The  present  procedure  is  ad- 
mittedly approximate,  since  it  is  essentially  based  on  a linearization  of  the 
nonlinear  interference  process,  but  it  appears  to  provide  an  adequate  repre- 
sentation of  the  flow-phenomena  (with  the  same  accuracy  of  more  complex  anal- 
yses) at  least  for  the  cases  reported  here. 

A shock  fitting  technique  has,  in  addition,  been  developed  for  the  exact  deter- 
mination of  shock  location  and  shape  over  a supercritical  airfoil.  A concomit- 
ant computer  program  has  been  developed  and  has  been  used  to  precisely  assess 
the  effect  of  a perturbation  of  back  pressure  on  shock  location. 
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APPFIJDIX  A 

BLOCKAGE  INTERFERENCE  BY  LINEARI7ED  SUBSONIC  THEORY 


The  linearized  subsonic  theory  of  wall  interference  for  nonlifting,  tv/o- 
dimensional  airfoils  in  tunnels  with  solid,  porous  and  slotted  walls  is  re- 
viewed here.  The  present  derivation  determines  the  interference  potential  by 
following  the  analysis  of  Lo^'^^  for  the  particular  case  where  the  non-linear 
term  is  neglected. 


The  analysis  is  based  on  the  solution  of  the  linearized  equation  of  subsonic 
compiess i bl e flow 


(1 


XX 


yy 


= 0 


(A.l) 


subject  to  the  tangential  flow  condition  at  the  airfoil  surface 


v;hcre  F{x)  is  the  airfoil  profile  function,  and,  for  the  case  of  a porous  wall, 
to  the  wall  boundary  condition,  at  y = h, 

❖x  + T <t,y  = 0 (A. 3) 

where  T is  the  porosity  parameter.  For  T ->■  <»  the  solid  v/all  boundary  condi- 
tion, = 0,  is  recovered. 


The  method  of  Fourier  transform  is  now  applied,  by  taking  the  Fourier  trans- 
form of  ❖ (x) , 


(|i(p)  = (27i) 
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(j>(x)  dx, 
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to  yield  the  equation 
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and  ttie  boundary  conditions 


= T f(p)  = (2,,)' 


T e dx,  at  y = 0 


- ip  <!)  + T (j,  = 0, 

y 


at  y = 0 


The  solution  of  Equation  (A. 5),  subject  to  boundary  conditions  (A. 6)  and  (A.?), 
in  the  transform  plane  is 

-Xh 

^ = A - LJApI  (ip  + TX)e 

n X ip  cosh  Xh  + tx  sinh  Xh 

where  X = Pp  and  where  is  the  free  air  solution  i.e.,  = - (t  f(p)/X)e’^^. 

The  interference  potential  in  the  transform  plane,  is  then  given  by 


({;  - 4.  - (} 

* m 


which,  fy  the  inverse  transform,  yields,  in  the  physical  plane,  the  real  value 
potential 


“Xh  cosh  X\ 


I “Xh  cosh  Xy 

dp  I e lb  M 

s AS 
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where 

1 

c 

= (p^  H + TX  p^G)/D 

's 

= (-TX  M + P^G)  p/D 

D 

2 2 2 
= H + p'^  G 

H 

= TX  sinh  Xh 

G 

= cosh  Xh 

X = fSp 


and  where  and  are  related  to  the  airfoil  geometry  through  the  integrals 


M 
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^ F(f,)  cos  p (x-4)  d^ 
-1/2 


(A. 9) 


and 
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F(5)  sin  p (x-fj 


'-1/2 


(A. 10) 


In  the  case  of  a tunnel  with  slotted  walls,  the  wall  boundary  condition  is 
usually  given  by  the  homogeneous  form 


((1  + K 41  = 0 

X xy 


(A. 1 1 ) 


where  K is  the  slot  parameter.  For  convenience,  the  integrated  form  of  Equa- 
tion (A. 11)  is  used  here,  i.e.,  the  form 
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(|)  + K(f^  = 0 aty  = h 

y 

which  implies  zero  perturbation  potential  at  upstream  infinity.  Again  by  us- 
ing the  Fourier  transform  (A.'i),  the  solution  of  Equation  (A. 5)  subject  to  the 
boundary  condition  (A. 6)  at  the  airfoil  surface,  and  to  the  boundary  condition 


({«  + K 4'y  = 0 at  y = h 


yields,  in  the  transform  plane,  the  interference  potential 


_ j (l-X  K)e  2 cosh 


(HX  K)c^'’  + (1-X  K)e"^*’ 


By  the  inverse  transfrom,  we  then  obtain,  in  the  physical  plane,  the  interfer- 
ence potential 


where 


T I.  „(i“XK)e^^2cosfiXy 

II  I*  ^ Xh  M dp 

' (1+X  K)c^  + (1-X  K)e 
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